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ABSTRACT 

Wc report the discovery of a large-amplitude oscillation in the hot subdwarf B star CS 
1246 and present multi-colour photometry and time- resolved spectroscopy supporting 
this discovery. We used the 0.41-m PROMPT array to acquire data in the u', g', r', and 
i' filters simultaneously over 3 consecutive nights in 2009 April. These data reveal a 
single oscillation mode with a period of 371.707 ± 0.002 s and an amplitude dependent 
upon wavelength, with a value of 34.5 ± 1.6 mma in the u' filter. We detected no 
additional frequencies in any of the light curves. Subsequently, we used the 4.1-m 
SOAR telescope to acquire a time-series of 248 low-resolution spectra spanning 6 hrs 
to look for line profile variations. Models fits to the spectra give mean atmospheric 
values of T^ef f = 28450 ± 700 K and log g — 5.46 ± 0.11 undergoing variations with 
semi-amplitudes of 507 ± 55 K and 0.034 ± 0.009, respectively. We also detect a 
radial velocity oscillation with an amplitude of 8.8 ± 1.1 km s^^. The relationship 
between the angular and physical radii variations shows the oscillation is consistent 
with a radial mode. Under the assumption of a radial pulsation, we compute the stellar 



distance, radius, and mass as d 



460 ± 11° 



pc, R = 0.19 ± 0.08 R0, and M = 0.39 



±0 ^3 Mq, respectively, using the Baade-Wesselink method. 
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1 INTRODUCTION 

Hot subdwarf B (sdB) stars are evolved objects presumed to 
have He-burning cores and modelled as extended horizontal 
branch stars (Heber 1986). They dominate surveys of faint 
blue objects and are often cited as the main source of the 
UV-upturn observed in globular clusters and giant elliptical 
galaxies. Their temperatures range from 20000 to 40000 K, 
and their log g values span more than an order of magnitude 
from 5 to 6.2. Their evolutionary histories are not entirely 
understood; various formation scenarios are possible includ- 
ing both single-star and binary-star channels (D'Cruz et al. 
1996; Han et al. 2002, 2003). At least two thirds of them 
appear to be in binary systems (Saffer et al. 2001; Maxted 
et al. 2001), but a higher binary fraction cannot be ruled 
out. Han et al. (2002, 2003) have even shown it is possible 
to explain 100% of the sdB stars as the products of binary- 
star channels. In either case, sdB stars probably evolve from 
the red giant branch after losing a signficant portion of their 
outer H envelopes, and stellar models show they will evolve 
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directly to the white dwarf cooling sequence after exhaust- 
ing the He in their cores (Dorman et al. 1993). Their optical 
spectra are dominated by H Balmer lines and sometimes also 
display He I or even the 4686 A He II line. See Heber (2009) 
for a review on what is currently known about sdB stars. 

More than a decade ago, Kilkenny et al. (1997) re- 
ported the first detection of photometric oscillations in a sdB 
star (EC 14026-2647). Contemporaneous with this discov- 
ery, Charpinet et al. (1996, 1997) predicted the existence of 
a class of sdB pulsators after calculating stellar models that 
showed driving of p-modes. The oscillations they found were 
driven by the K-mechanism and associated with an opacity 
bump from the ionisation of Fe in their model envelopes. 
The discovery of the first pulsator sparked an interest in 
sdB variability studies, and, subsequently, around 70 pul- 
sators have been discovered to date. Many techniques have 
been applied to sdB pulsators successfully including multi- 
colour photometry (see Tremblay et al. 2006 for a summary) 
and time-resolved spectroscopy (pioneered by O'Toole et al. 
2000; Telling & 0stensen 2004; TiUich et al. 2007). Astero- 
seismological studies of these stars have revealed information 
about their internal structure, including total mass, envelope 
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mass, and chemical stratification and have been carried out 
on a handful of sdB pulsators (see Brassard et al. 2001 and 
Charpinet et al. 2005 as examples). 0stensen (2009) presents 
a thorough review of asteroseismolgy and its successes. 

Obscrvationally, the subdwarf B pulsators may bo di- 
vided into three classes: the rapid (sdBVr), slow (sdBVs), 
and hybrid (sdBVrs) pulsators^. The sdBVr stars are as- 
sociated with models of p-mode pulsators and have periods 
between 80 and 600 s and typical amplitudes less than 1%. 
In the log g-Teff plane, they cluster together in an insta- 
bility strip with Te// ranging from 28000 to 35000 K and 
log g between 5.2 and 6.1. The sdBVs stars generally have 
periods from 1 to 2 hrs with amplitudes of only a few tenths 
of a percent; their oscillations are modelled with g modes. 
They are also found in an instability strip with Teff be- 
tween 23000 and 30000 K and log g near 5.4. In the region 
where the blue edge of the sdBVs instability strip overlaps 
the red edge of the sdBVr strip, three sdBVjs stars have 
been observed with fast and slow oscillations (Schuh et al. 
2006; Baran et al. 2006; Lutz et al. 2009). Both of the in- 
stability strips contain stars that have been observed not to 
vary photometrically. 

The most fundamental parameters of any star, the mass 
and radius, are not well known for sdB stars and are rarely 
computed from observations. Single and binary evolution 
calculations point to a canonical mass of 0.47 Mq (Han et al. 
2002, 2003), and a few masses supporting this result have 
been derived from astcroscismology (e.g.. Brassard ct al. 
2001; Randall et al. 2007) and from eclipsing binaries (e.g.. 
Wood et al. 1993, Drechsel et al. 2001). Zhang et al. (2009) 
have even proposed a new technique for deriving statistical 
masses from observational data in the literature. Most of the 
sdB radii measurements to date come from binary systems 
for which the separation distance is known; the sdB radius 
is derived from the ratio of the radii and assumes informa- 
tion about the companion. Such studies show a typical sdB 
radius to be near 0.2 Rq, but additional and more direct 
measurements are needed to confirm this result. 

We report the detection of a single, laxge-amplitude os- 
cillation in the sdB star CS 1246 (q;2000=12:49:37.69, (520oo = 
-63:32:08.99), an object first observed by Reid et al. (1988) 
in a survey of hot white dwarfs and OB subdwarfs in ob- 
scured regions of the Galactic plane. In one of the most 
opaque sections of the Coalsack Dark Nebula, CS 1246 was 
identified by Reid et al. (1988) as a candidate hot subdwarf 
B star with a V magnitude of 14.59 ± 0.14 and U-B and 
B-V colours of -1.08 ± 0.10 and 0.28 ± 0.10, respectively 
Neither T^ff nor log g has previously been determined. 

We first learned of CS 1246 while browsing the on- 
line Subdwarf Database (0stensen 2006) for candidate pul- 
sators observable with the 0.41-m Panchromatic Robotic 
Optical Monitoring and Polarimetry Telescopes (PROMPT) 
in Chile. Even though there were no available estimates of 
Teff or log g, we decided to look for pulsations since the 
U-B colour roughly matched that of known pulsators and 

^ Here wc adopt the naming scheme in which all sdB pulsators 
are given the base name 'sdBV; subscripts 'r' and 's' arc added 
upon observations of rapid and slow oscillations, respectively. For 
reference, the sdBVr stars are also called V361 Hya and EC 14026 
stars, while the sdBVs stars have the aliases PG 1716, V1093 Her, 
and 'Betsy' stars. 



the star was a brighter, Southern-hemisphere target. The 
discovery run data revealed an interesting light curve, and 
we subsequently obtained multi-colour time-series photome- 
try and time-resolved spectroscopy to study the pulsations. 
We observed line profile variations in the spectra and found 
substantial oscillations in the radial velocity, temperature, 
and gravity. The nature of these variations is consistent with 
that expected for an £ = pulsation, but we cannot rule out 
the possibility of non-radial modes with our current data 
set. Although not one of our original goals, wo were able to 
compute the distance, radius, and mass of CS 1246 from the 
relation between the physical and angular radii variations. 
Our study is the first to apply such a method to an sdBV 
star, and the results show the mass and radius to be con- 
sistent with those calculated for other sdB stars in binaries 
and from theory. 



2 DISCOVERY OF THE PULSATION 

We first observed CS 1246 on 2009 March 30 after sub- 
mitting a job to the PROMPT queue using the web-based 
SKYNET interface (see §3.1 for a description of these sys- 
tems). Later that night PROMPT 4 monitored the star for 
approximately 75 min, obtaining an uninterrupted series of 
80 exposures through a Johnson V filter with integration 
times of 40 s and a duty cycle of 89%. Even before removing 
atmospheric extinction effects or transparency variations, wc 
realized CS 1246 was a pulsator as its oscillations were visi- 
ble in the undivided, raw light curve. After formally reducing 
and analysing the data using the methods described in §3.2 
and §3.3, we found a single oscillation with a period of 372 
s and an amplitude of 23 mma. Figure 1 shows the reduced 
light curve and its amplitude spectrum. No other frequencies 
were observed in this data set, which had a mean noise level 
of 1.6 mma. The relatively long period and large amplitude 
of the oscillation inspired us to obtain simultaneous multi- 
colour photometry with the PROMPT and timc-rcsolvod 
spectroscopy with the SOAR telescope. These observations 
and the results thereof are described in the sections that 
follow. 

3 MULTI-COLOUR PHOTOMETRY 
3.1 The instrument 

PROMPT is a system of 0.41-m Ritchey-Chretien telescopes 
at Cerro Tololo Inter-American Observatory (CTIO) in 
Chile built primarily for rapid and simultaneous multiwave- 
lengtli observations of gamma-ray burst afterglows. Table 1 
presents some of the CCD characteristics of these telescopes. 
The PROMPT telescopes are 100% automated and operate 
under the control of SKYNET, a prioritised queue schedul- 
ing system that allows registered users to acquire, view, and 
download data from any location via a PHP-enabled Web 
server. Observations in the queue are assigned priority lev- 
els and scheduled accordingly; gammarray bursts receive top 
priority. The SKYNET system enables automated data col- 
lection at times that would be inconvenient or impossible for 
a human observer. For example, our multi-colour photomet- 
ric programme was executed by PROMPT while BNB and 
BHD were in Detroit, Michigan, supporting the University of 
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Figure 1. Light curve (upper panel) and amplitude spectrum 
(lower panel) of the discovery run data from 2009 March 30. The 
light curve was produced from V-filtered data and is presented 
unsmoothed. 



Table 1. Characteristics of the PROMPT CCDs 



Telescope 


Size 


Pixel Scale 


Gain 


Rcadnoisc 




(pix) 


(arcsec/pix) 


(e-/ADU) 


(e-) 


PROMPT 2 


3k X 2k 


0.41 


1.6 


8.3 


PROMPT 3 


Ik X Ik 


0.59 


1.6 


8.3 


PROMPT 4 


Ik X Ik 


0.59 


1.4 


7.6 


PROMPT 5 


Ik X Ik 


0.59 


1.4 


9.2 



North Carolina Men's Basketball team at the 2009 NCAA 
National Championship. Obtaining the data presented in 
§3.2 using classical observing methods would have been an 
expensive and time-consuming undertaking. For additional 
details on PROMPT and SKYNET, refer to Reichart et al. 
(2005). 



3.2 Observations and reductions 

PROMPT obtained multi-colour photometry simultane- 
ously through the u', g', r', and i' filters on 2009 April 4, 
5, and 6 using all available telescopes. Table 2 presents a 
detailed log of these observations. The CCD parameters of 
the PROMPT telescopes (Table 1) and SKYNET overhead 
resulted in duty cycles of 94%, 80%, 89%, and 89% for the 
u', g', r', and i' filter data sets, respectively, during the con- 
tinuous portions of the runs. Although the majority of the 
runs were uninterrupted, each data set contains an approx- 
imately 15-minuto long gap duo to the inability of the cq- 
uitorial mounts to track through the meridian. Additional 
40-minute gaps are present in the u' and g' data from 2009 
April 6. In total, our 3-day photometry campaign produced 
113 hours of data with 7291 useable frames. 

After bias-subtracting, dark-subtracting, and flat- 
fielding the frames with IRAF using conventional proce- 
dures, we extracted our photometry using the external IRAF 
package CCD_HSP written by Antonio Kanaan. Aperture 
radii were chosen to maximize the signal-to- noise (S/N) ra- 
tio in the light curves and were always less than the seeing 
width. We used sky annuli to subtract the sky counts from 
each stellar aperture, being sure to avoid contamination of 
the annuli with other stars in the field. Wo divided the sky- 
subtracted light curves by an average of constant comparison 
star light curves to correct for transparency variations; these 
stars are labeled in the field image shown in Figure 2. Some 
of them lacked an adequate S/N ratio through some filters 
while saturating the CCD in others, and, consequently, we 
employed different comparison stars for different filters (see 
Table 2). We are confident that the results of our frequency 
analyses are independent of the choices of the comparison 
stars, as they all appeared to be constant. Atmospheric ex- 
tinction efi'ects were approximated and removed each night 
by fitting and normalising the curves with parabolas. We 
used the WQED suite (Thompson & Mullally 2009) to ana- 
lyze the data and apply timing corrections to the barycentre 
of the solar system. Figure 3 presents the full, differentially- 
corrected light curves. 
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Figure 2. CS 1246 field. Image shown is a stack of 30 40-s V-band 
frames taken with PROMPT 4. CS 1246 is indicated with a star 
while the comparison stars are marked with circles and labeled 
C1-C9. 



3.3 Frequency analysis 

We analysed our reduced light curves using two tools: the 
discrete Fourier transform and the least-squares fitting of 
sine waves. WQED, again, was used to access these tools. 
Figure 4a presents the amplitude spectra for the combined 
light curves shown in Figure 3 along with their window func- 
tions. The u' and g' amplitude spectra are plotted out to 
their respective Nyquist frequencies, while the r' and i' spec- 
tra are terminated at 10 mHz (beyond this point these spec- 
tra are consistent with noise). 

Inspection of the amplitude spectra reveals a single, 
dominant mode at 2690.29 jjJiz (371.707s) with amplitudes 
of 32.0 ± 1.5, 24.4 ± 0.9, 21.3 ± 0.3, and 20.6 ± 0.4 mma 
in the u', g', r', and i' passbands, respectively. Table 3 sum- 
marizes the results of the least-squares fits; the errors shown 
are derived directly from the fits and assume uncorrelated 
noise. The results of Montgomery & O'Donoghue (1999) im- 
ply these errors may be larger by a factor of three. We detect 
no phase differences between the four colour curves within 
our level of precision. As the integration times of our expo- 
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Table 2. Multicolour photometry observations log 





UciLe 


T:^iif 


Telescope 


T 

-■- exp 


T 

cycle 


Length 


^ Frames 


Coinpa,rison Stars 


{V 1 J 


f'TTT^ 
1 J 










(hrs) 






23:59:02.8 


2009 Apr 04 


u' 


PROMPT 3 


80 


86 


9.4 


393 


C1,C2,C9 


00:57:03.3 


2009 Apr 05 


g' 


PROMPT 2 


40 


52 


8.7 


592 


C1,C2,C3,C4,C8,C9 


00:01:47.7 


2009 Apr 05 


r' 


PROMPT 4 


40 


47 


9.6 


710 


C2,C3,C4,C5,C7 


00:01:51.7 


2009 Apr 05 


i' 


PROMPT 5 


40 


47 


9.6 


710 


C2,C3,C4,C5,C6 


23:57:56.0 


2009 Apr 05 


u' 


PROMPT 3 


80 


86 


9.4 


380 


C1,C2,C9 


23:57:05.9 


2009 Apr 05 


g' 


PROMPT 2 


40 


52 


8.6 


580 


C1,C2,C3,C4,C8,C9 


23:57:41.8 


2009 Apr 05 


r' 


PROMPT 4 


40 


47 


9.3 


700 


C2,C3,C4,C5,C7 


23:57:38.9 


2009 Apr 05 


i' 


PROMPT 5 


40 


47 


9.3 


700 


C2,C3,C4,C5,C6 


23:56:41.8 


2009 Apr 06 


u' 


PROMPT 3 


80 


86 


9.7 


395 


C1,C2,C9 


23:55:47.2 


2009 Apr 06 


g' 


PROMPT 2 


40 


52 


9.7 


663 


C1,C2,C3,C4,C8,C9 


23:56:36.1 


2009 Apr 06 


r' 


PROMPT 4 


40 


47 


9.7 


735 


C2,C3,C4,C5,C7 


23:56:31.8 


2009 Apr 06 


i' 


PROMPT 5 


40 


47 


9.7 


733 


C2,C3,C4,C5,C6 



"Time at midpoint of first exposure. 
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Figure 3. Light curves of CS 1246 over 3 consecutive nights (left panels) taken simultaneously through the u', g', r', and i' filters (from 
top to bottom). These observations represent a total of 113 hours of data taken with 7291 frames. Expanded views (right panels) of the 
same 1.5-hour section of the curves are shown to illustrate the variable nature of the star. 



sures are significant fractions of tiie 371. 7-s period, we must 
talis pliase smearing into account. For any measured observ- 
able oscillating witii a period P and a physical amplitude A, 
a finite integration time / will result in a reduced observed 
amplitude Ao given by 

A, = —stn{—)A, (1) 

Baldry (1999) presents a discussion of phase smearing and a 
derivation of this equation. The amplitudes we measure for 
the 371. 7-s variation underestimate the physical amplitudes 
by 7.4% in the u' filter and 1.9% in the g', r', and i' fihers. 
The smearing-corrected amplitudes are 34.5 ± 1.6, 24.8 ± 



0.9, 21.7 ± 0.3, and 21.0 ± 0.4 mma in the u', g', r', and i' 
filters, respectively. 

After fitting the dominant mode and subtracting the fit 
from the light curves, we recalculated the Fourier transforms 
of the residual light curves; they are shown in Figure 4b. The 
mean noise levels in these pre-whitened u', g', r', and i' am- 
plitude spectra are 1.9, 1.1, 0.3, and 0.5 mma, respectively. 
All of the spectra show a heightened noise level at extremely 
low frequencies, but this is likely a side effect of our method 
for extinction removal. Beyond the low-frequency regime, 
some of the pre-whitened spectra show peaks at or above the 
3-cr level. However, none of these periodicities were detected 
in more than one data set, whether from night-to-night or 
from filter-to-filter, so we evaluate their signficance using 



Photometry and spectroscopy of CS 124 6 5 

(a) 



0.03 

0.02 

0.01 
0.00 

^ 0.02 

|- 0.01 

- 0.00 
§ 0.02 



2 0.01 - 



0.00 
0.02 



0.01 - 
O.OOQftttuu 



^ 1 



j»^». |. M | I ' 1 1. 1 



|,i| r j - I ,| „ | .. . | . | , I .| n |. ,.| ,), M |, I il | 



1 1 1 1 1 1 1 1 1 1 1 




i f-filVfff'v 1 1 1 


Mill 'I'-AVrffr 


: J 













0.010 — 



2000 



4000 6000 
Frequency {fiHz) 



8000 



-200 -100 100 200 
Frequency (nHz) 



(b) 



0.005 

0.000 
03 0.006 

T3 

i. 0.003 
E 

f 0.000 1 

CO 

o 




it' ■ 




4000 6000 
Frequency (|xHz) 



10000 



Figure 4. Amplitude spectra of tlie combined liglit curves in the u', g', r', and i' filters (from top to bottom) before (a) and after (b) 
pre- whitening of the main mode. The window functions for each data set are shown in the right portion of (a). In the pre- whitened 
spectra, the So" and Aa levels are represented by dashed and dot-dashed lines, respectively, and the position of the pre- whitened signal 
is marked with a vertical dashed line. The u' and g' spectra are plotted out to their respective Nyquist frequencies. The r' and i' plots 
are terminated at 10 mHz but are consistent with noise beyond this point. 



the false alarm probability (Horne & Baliunas 1986). None 
of them meet our significance test of a false alarm probabil- 
ity <1%. Moreover, the folded pulse shapes agree with that 
of a single sinusoid to within the errors, implying a lack of 
harmonically-related frequencies in the light curve. 

The amplitude of the oscillation is strongly wavelength- 
dependent, as illustrated in Figure 5. The u' amplitude is the 
largest, being 39%, 59%, and 64% greater than those in the 
g', r', and i' filters, respectively. In principle, these ratios may 
be used along with the known atmospheric parameters of the 
star to estimate the degree / of the oscillation, but in our 
case the large error bars make such an exercise inconclusive. 



4 TIME-RESOLVED SPECTROSCOPY 
4.1 The instrument 

The Goodman Spectrograph is an imaging spectrograph 
built for one of the Nasmyth ports of the 4.1-m SOAR 
telescope on Cerro Pachon in Chile. In imaging mode, the 
camera-coUimator combination re-images the SOAR focal 
plane with a focal reduction of three times, yielding a plate 
scale of 0.15 arcsec pixel" ^ at the CCD. The camera con- 
tains a 4k X 4k Fairchild 486 back-illuminated CCD with 
electronics and dewar provided by Spectral Instruments, Inc. 
of Tucson, Arizona. The system uses optics of fused silica. 
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Table 3. Results of the least-squares fits to the light curves 



Filter 


Frequency 


Amplitude'' 


Phase= 




(mHz) 


(mma) 


(s) 


u' 


2690.28 ± 0.10 


32.0 ± 1.5 


76.5 ± 2.8 


g' 


2690.30 ± 0.08 


24.4 ± 0.9 


75.7 ± 2.2 


r' 


2690.29 ± 0.03 


21.3 ± 0.3 


77.3 ± 0.8 


i' 


2690.30 ± 0.05 


20.6 ± 0.4 


77.5 ± 1.2 



^ Not corrected for phasc-smcaring. 

" Time of maximum measured from BJED 2454925.50315492. 




3000 4000 5000 6000 7000 8000 
Wavelength (A) 

Figure 5. Amplitude ratios A^//A^/ of the 371.707s mode for the 
u', g', r', and i' data, corrected for phase-smearing. Normalisation 
to A^i was chosen to minimize the error bars. 



NaCl, and CaF2 to achieve high throughput from 320 to 850 
nm. In spectroscopic mode, the instrument can obtain spec- 
tra of objects over a 3'x 5'field using long slits with widths 
ranging from 0.46" to 10". A multi-slit mode is currently 
being commissioned. Users may choose from three Volume 
Phase Holographic (VPH) gratings with dispersions of 0.31, 
0.65, and 1.3 A pixel For further details concerning the 
spectrograph, refer to Clemens et al. (2004). 



4.2 Observations and Reductions 

We obtained 248 low-resolution spectra of CS 1246 during an 
engineering run on 2009 April 17 with the Goodman Spec- 
trograph in order to look for line profile variations. The 600 
mm~^ VPH grating (0.65 A pixol"^ dispersion) was used to 
cover a spectral range from 3500 to 6200 A. As the sky at 
Cerro Pachon was photometric with reasonable seeing, we 
employed a 10"-wide slit and let the average seeing of 1.1" 
determine our final resolution of 4.8 A. The position angle 
was set to 272.3 degrees E of N so that we could place a 
nearby comparison star on the slit. We binned the spectral 
images by two in the spatial and dispersion directions and 
only read out 400 binned pixels along the slit to minimize 
the readout time. Each exposure had an integration time 
of 80 s with 6 s of overhead between images, resulting in 
a duty cycle of 93%. The entire run lasted 6 hours from 
00:59:42.0 to 07:07:53.4 UT with a 15-minute gap that be- 
gan near 02:19:00 UT. Comparison spectra of HgAr and Cu 
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Figure 6. Individual (top) and average (bottom) spectra of CS 
1246 taJsen with the Goodman spectrograph on 2009 Apr 17. The 
spectra shown were flux-calibrated using the spectrophotometric 
standard EG 274. 

lamps were obtained for wavelength calibration during the 
gap and after the entire series was complete. 

We bias-subtracted and flat-ficldcd our spectra using 
standard tasks within IRAF. The apall routine was used to 
optimally extract one- dimensional spectra from the frames 
and subtract a fit to the sky background. The resulting in- 
dividual spectra have an average S/N ratio of 49 pixel" ^ 
in the continuum near 5000 A. The top panel of Figure 6 
shows an example of one of these spectra. Unfortunately, we 
were not succesfuUy able to use the HgAr or Cu frames to 
wavelength-calibrate our spectra, as most of the 43-A wide 
lines overlapped and made it difficult to accurately mark line 
centres. As an alternative, we self-calibrated the mean spec- 
trum to the H Balmer and He I lines and used the solution to 
wavelength-calibrate the individual spectra. Although this 
technique makes it impossible to measure the true space ve- 
locity of CS 1246, it preserves the relative wavelength scale 
without removing offsets between spectra and does not efi^ect 
our relative radial velocity measurements. We checked the 
resulting wavelength solution with that derived from past 
spectra taken with the same instrumental setup (but with 
smaller slits) and found that the dispersions agreed well. 

4.3 The mean spectrum 

The average spectrum has a S/N ratio of 770 pixel" ^ and is 
shown in the bottom half of Figure 6. H Balmer (H/3-H13) 
and He I lines are the most prominent features in the spec- 
trum, but signatures of Si HI, N II, and O II are also present. 
The spectrum appears to be contaminated with Ca II (3933 
A) and Na D. Although these features usually indicate the 
presence of a cooler companion, we expect some level of 
interstellar absorption given the extremely low galactic lat- 
itude (-0.94 deg) of CS 1246. Moreover, both Ca and Na 
D have been observed in the spectra of other stars located 
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in and behind the Coalsack (Franco 2000, 1989; Crawford 
1991). No other spectral features indicative of a companion 
are present, and we currently have no reason to believe the 
lines originate from a companion star. 

To determine the atmospheric parameters from the 
mean spectrum, we employed the model grid of synthetic 
optical spectra described in Han et al. (2007). This grid is 
defined in terms of 21 values of T^/ / from 20000K to 40000K 
in increments of lOOOK, 11 values for the surface gravity 
from log g of 5.0 to 7.0 in increments of 0.2, and 4 values 
of the helium-to-hydrogen ratio from log Ar(He)/iV(H) of - 
3 to in increments of 1. To obtain further resolution in 
the grid, we interpolated the models using a cubic spline 
interpolation technique (see Press et al. 2007). Each model 
spectrum was degraded by convolving it with a Gaussian 
whose FWHM mimicked that of a single resolution element 
in our spectrum. Finally, we re-binned the model spectra to 
match the binning of the observed spectrum. 

The comparison to the model spectra needed to be car- 
ried out in terms of normalised line profiles, and so we im- 
plemented a fitting technique similar to that described in 
Bergeron et al. (1995). We normalised the flux in the He 1 
and H Balmer (H/3-H9) profiles to a linear continuum set to 
unity at two points. These points are located at fixed dis- 
tances on either side of the profile and are far enough from 
the line centre to completely encompass the absorption pro- 
file. The unity point fluxes can be determined by averaging 
the fluxes in the neighbouring pixels, but defects in the ob- 
served spectrum at these locations may lead to inaccurate 
fits to the continuum. Instead, we fit each absorption line 
and its surrounding continuum with a "pseudo" Lorentzian 
profile as described in Saffer et al. (1988). This function has 
the form 
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Figure 7. Light curve (top) produced from the spectrophotom- 
etry in the bandpass 3500-6200 A. Although a paraboUc fit has 
been removed to take out extinction variations, the curve has not 
been divided by that of a constant comparison star, attesting to 
the superb photometric quality of the night. The amplitude spec- 
trum (bottom) of the curve has a mean noise level of 0.65 mma. 



The exponent a controls the overall shape of the function 
and fits both the peak and the wing of an absorption line 
much better than a classic Lorenztian or Gaussian profile. 
We used the IDL routine MPFIT (Markwardt 2009), which 
employs the Levenberg-Marquardt method (Press et al. 
2007), to perform a non-linear least-squares fit to each line. 
The unity point fluxes are derived from the resulting smooth 
fit, and each profile is normalised by dividing by the straight 
line passing through these points. We normalised the line 
profiles in model spectra in the same manner except that 
the unity point fluxes were taten directly from the spec- 
tra (which are noiseless) instead of being determined by a 
smooth fit. 

We determined the best-fit model for each line profile 
with a minimization technique and found Tef / = 28450 
± 300, log g = 5.46 ± 0.06, and log N{Rc)/N{}l) = -2.0 
±0.1. Systematic errors in the data and from fitting to im- 
perfect models also contribute to the uncertainties in these 
values. Consequently, we believe the statistical errors pre- 
sented above underestimate the uncertainties and adopt the 
following values: Teff = 28450 ± 700, log g = 5.46 ± 0.11, 
and log iV(He)/Af(H) = -2.0 ± 0.3. These parameters place 
CS 1246 near the boundary of the sdBVr- and sdBVs instar 
bility strips. 



4.4 Spectrophotometric light curve 

We created a light curve from our spectroscopic data by 
summing together the flux from all wavelength bins in each 
of the 248 spectra. To correct for transparency variations, 
we attempted to divide the resulting light curve by that of 
the comparison star, but its low S/N ratio added noise to the 
light curve, and we abandoned the idea. We did, however, 
correct for atmospheric extinction effects by fitting and nor- 
malising the light curve with a parabola. We note that this 
normalisation removes any intrinsic variations on the order 
of the run length. In spite of the absence of a comparison 
star by which to divide, the high photometric quality of the 
night and good S/N ratio in the spectra led to the superb 
light curve shown in Figure 7 above its FT. 

A least-squares fit to the curve reveals a single frequency 
of 2689.9 ± 0.6 ^iHz (371.76 ± 0.09 s) with an amplitude of 
21.5 ± 0.6 mma. As with the time-series photometry, we 
correct for phase smearing, which leads to a damping of 
7.45%, and derive an amplitude of 23.3 ± 0.7 mma. We de- 
tect no additional periodicities with signficant probabilities 
from the light curve, which had a mean noise level of 0.65 
mma. To compare the spectrophotometry to the PROMPT 
results, we convolved the spectra with the g' filter trans- 
mission function and repeated the light curve analysis. We 
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Figure 8. Radial velocities obtained from fitting pseudo 
Lorentzian profiles to the stellar H/3 line (top) and the 5577 A 
O I sky emission line (bottom). The large-scale variation is due 
to instrumental flexure. 



find a smearing-corrected amplitude of 23.2 ± 0.7 mma in 
agreement with the g' filter result from §3.3 . 



4.5 Radial velocity curve 

To construct a radial velocity curve, we individually fit the 
H/3-H9 profiles in each spectrum with "pseudo" Lorentzian 
functions as described in §4.3, fixing the a parameter to the 
value determined from the average spectrum. The individual 
spectra did not have the S/N ratio needed to fit the other 
absorption line profiles accurately, and, consequently, these 
lines were ignored in our analysis. The locations of the profile 
centres in the mean spectrum were subtracted from those in 
the individual spectra, and the shifts from the mean were 
converted to km s~^. 

The top panel in Figure 8 shows the resulting radial ve- 
locity curve for H/3, as an example. The data exhibit a long- 
term variation spanning approximately 100 km s^^ over the 
course of the run. We attribute this structure to the com- 
bination of two factors: the global change in instrumental 
fiexure and the movement of the star at the focal plane. As 
the sky fills the entire slit, we may correct the former by ex- 
amining the drift of atmospheric emission lines. The latter 
problem is more difficult to correct, and, consequently, we 
do not attempt to do so. Previous time-series photometry 
runs with the Goodman Spectrograph, however, show the 
drift due to guiding to be small, usually around 0.3" over a 
time span much longer than the period of CS 1246. 

To compensate for flexure, we flt the O I sky line in each 
spectrum with a Gaussian using the GAUSSFIT routine in 
IDL and obtained the curve shown in the bottom panel of 
Figure 8. A 5th-order polynomial matched the O I curve 
fairly well, and we removed this trend from the velocity 
curves of the stellar absorption lines. After fiexure correc- 
tion, the velocity curves were linear with slopes dependent 
upon wavelength. This result is expected from atmospheric 
refraction and allows us to approximate the guiding wave- 
length as 8400 A. We computed the global radial velocity 
curve by removing the linear variations due to refraction 
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Figure 9. Flexure-corrected radial velocity curve (top) produced 
from the weighted average of the H/3-H9 velocity curves. The am- 
plitude spectrum (bottom) has a mean noise level of 0.94 km s^^. 
Shortcomings in our flexure-correction are visible in the 3.25-4.0 
hr section of the curve and give rise to the low-frequency peak in 
the FT near 143 /xHz. 



and averaging the velocities from all of the H Balmer lines, 
weighting each curve by the inverse of its variance. Figure 9 
shows the resulting radial velocity curve along with its FT. 

A strong radial velocity variation accompanies the pho- 
tometric oscillation at the same frequency. After barycen- 
tric correction of the velocities and the mid-exposure times, 
we find a semi-amplitude of 8.2 ± 1.0 km s~^ by fitting 
a sine curve with frequency fixed to that derived from the 
PROMPT data. Sine wave fits to the velocity curves of the 
individual H/?-H9 lines give amplitudes of 8.7, 7.2, 7.5, 8.6, 
8.9, and 9.8 km s~^, respectively. The uncertainty we re- 
port for the average velocity is the standard deviation of 
these values. The increase in velocity with increasing Balmer 
order (with the exception of H/3) implies the velocity is 
wavelength-dependent, as expected from weighting of the 
projected velocities over the observed hemisphere by limb- 
darkening. Correcting for phase smearing, we derive a final 
semi-amplitude of 8.8 ± 1.1 km s~^. We claim no additional 
candidate frequencies after pre-whitening the data with the 
main pulsation period. The large peak in the FT near 143 
ytiHz is a result of an imperfect flexure correction and should 
not be attributed to stellar oscillations. The mean noise level 
in the amplitude spectrum is 0.94 km s"'^. 
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4.6 Effective temperature and gravity variations 

Wc determined Te// and log g as a function of pulsation 
phase by phase-folding and averaging the spectra together 
in 10 evenly-spaced bins. The resulting spectra had S/N rar 
tios of approximately 240 pixel" ^ in the continuum near 
5000 A. Wc fit atmospheric models to the spectra in the 
maimer described in §4.3 to determine Te//, log g, and 
log A'^(lle)/JV(H) over a pulsation cycle. As wc observed no 
change in the He abundance over the 10 phase bins, we fixed 
the abundance to its average value and refit the other atmo- 
spheric parameters. 

Figure 10 shows the Te// and log g variations plotted 
above the radial velocity and light curves, which were phase- 
folded into the same 10 bins. Since the temperature and 
gravity fits displayed sinusoidal variations, we fit sine waves 
to these curves using a fitting routine; the frequencies 
of the fits were fixed to the pulsational period. We derive a 
temperature variation of 462 K about an equilibrium value 
of 28279 K and a gravity variation of 0.031 about a value of 
5.44. These amplitudes are subject to phase smearing from 
two sources. The finite integration time of the individual 
spectra in the series results in an underestimation of the am- 
plitude by 7.45%. We must also consider the smearing that 
occurs when phase-folding and re-binning the spectra. Each 
phased-binned spectrum has a pseudo-integration time of 
37.17 s and leads to a further damping of 1.64%. Combiiung 
both smearing factors, we find that our fits to the Te/ / and 
log g curves underestimate the true variations by 8.97%. Ap- 
plying these corrections, we report semi-amplitudes of AT 
= 507 ± 55 K and Alog g = 0.034 ± 0.009. The errors are 
derived using the deviations from the fits; they are more 
precise than those associated with the absolute temperature 
and gravity since systematic errors tend to cancel out in a 
differential analysis. 

The phase relationships among the photometric, radial 
velocity, gravity, and temperature variations axe apparent in 
Figure 10. If conditions are purely adiabatic, the tempera^ 
ture would be 90 deg out of phase with the radial velocity 
(180 deg out of phase with the radius). Wo measure a dif- 
ference of -72 ± 22 deg. Given the size of the error bars, 
we cannot claim deviations from adiabatic conditions us- 
ing our current data set, although we do note this value is 
fairly deviant from 90 deg. The temperature, gravity, and 
flux variations are in phase to within the errors. 



4.7 Radial pulsation test and the stellar distance 

The simplicity of the light curve of CS 1246 combined with 
our high S/N spectra gives us the unique opportunity to 
check whether our observations are consistent with a radial 
oscillation. The method we use is similar to that described 
by Lynas-Gray et al. (1984) and relies on the computation of 
the angular radius a. In §4.6 wc found the best-fitting model 
spectrum to each of our observed phase-binned spectra, giv- 
ing us the temperature T at each phase (p. We take the 
monochromatic flux near 5500 A in the observed and model 
spectra for our calculations, as the continuum in this region 
is fairly smooth and free of absorption lines. The model flux, 
Fx , represents the flux at the photosphere and is a function 
of T only, by definition. The observed flux, fx, is a function 
of the temperature, the radius R, and the distance d. Con- 
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Figure 10. Variations in tlic effective temperature, surface grav- 
ity, radial velocity, spectropliotometry, and angular radius (from 
top to bottom) for the 371. 7-s mode. Non-linear, least-squares fits 
to the curves are shown as dotted lines, and the pulsation cycle 
is presented twice for visual clarity. The curves shown are taken 
directly from the data and have not been corrected for phase- 
smearing or projection effects. 



sequently, we can derive the angular radius a for each phase 
bin by computing the ratio of the model flux to the observed 
flux via the expression 



Fx{T) 



(3) 



fx{T,R,d) 

Flux calibration errors, atmospheric extinction, and inter- 
stellar extinction will affect the observed flux; the latter 
two lead to an underestimation of fx- We have already flux- 
calibratcd the spectra and corrected them for atmospheric 
extinction in our reductions, but we have not yet considered 
the effects of interstellar extinction. The temperature de- 
rived from our time-averaged spectrum implies a B-V colour 
of -0.29 ± 0.01. Reid et al. (1988) reported a B-V colour 
of -1-0.28 ± 0.10, from which we compute the reddening as 
E(B-V) — 0.57 ± 0.11. Assuming there is no companion star 
contributing to the reddening, we convert the E(B-V) esti- 
mate to an approximate extinction level of 1.8 mag for 5500 
A using the results of Schlcgcl et al. (1998). Thus, our ob- 
served fluxes should be larger by a factor of 5.2 (=2.5^ *). 
Taking this factor into account, we computed the angular 
radius for each phase bin and derived the a variation shown 
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in the bottom panel of Figure 10. A least-squares fit to the 
curve gives a semi- amplitude of 4.88*10"" rad (5.36*10"" 
rad after smearing-correction) about a mean angular radius 
of 1.01 * lO"'^'^ rad. We stress again that errors in the flux 
calibration and extinction estimates lead to incorrect values 
for the angular radii and their relative differences. 

A plot of Aa vs AR will reveal a linear relationship for 
a radial pulsation mode, as shown by the relation 

Aa=—, (4) 

where the distance is the inverse of the slope. Having already 
derived the a variations, we compute the variations in the 
physical radius by integrating the radial velocity curve. We 
converted the observed expansion velocities into actual ex- 
pansion velocities using a projection factor dependent upon 
relative limb darkening in the spectral lines. Montanes Ro- 
driguez & Jeffery (2001) computed such factors for early- 
type, radially pulsating stars, and their studies show a value 
of 1.4 to be appropriate for an sdB star like CS 1246. By 
applying this factor, we assume from this point on that the 
371. 7-s mode is radial in nature and can only check whether 
our observations are consistent with this assumption. The 
velocity amplitude becomes 12.3 ± 1.4 km s"^ after applica- 
tion of the projection factor, and integration of this velocity 
gives a physical radius variation with semi-amplitude of 728 
± 83 km. 

Figure 11 presents a plot of Aa against AR. To test for 
a linear relationship, we computed the Pearson correlation 

coefficient and find a value of 0.78. The probability that 10 
measurements of two uncorrelated variables can produce a 
correlation coefficient at or above this value is less than 2% 
(see Appendix A). We can neither rule out the possibility 
that the 371. 7-s mode in CS 1246 is a radial one, nor can 
we exclude the possibility of a mode with £ > 0. 

Still under the assumption of a radial mode, we fit a line 
to the data and determined the distance from the slope via 
Eqn. 4. As Aa = should correspond to the point where AR 
— 0, we force the fit to pass through this point and derive a 
distance of d = 460 ± loo Pc An sdB star this close should be 
brighter than what Reid et al. (1988) measured. In fact, with 
Te//=28450±700 K and R=0.19±0.08 Rq (see §4.8), CS 
1246 should have an apparent Visual magnitude of 12.7±[];5 
when observed at the above distance. This value implies a 
Visual extinction level of 1.9±o;e mags, in agreement with 
that found by Rodgers (1960) and Mattila (1970) for the 
part of the Coalsack where CS 1246 resides. 



4.8 The stellar radius and mass 

We are now in the position to compute the radius and mass 
directly from our observations using two different methods. 
The observed variations in log g come about from the pul- 
sational acceleration and changes in the radius, and after 
the former contribution is subtracted, the radius may be 
derived from the remaining difference. Unfortunately, the 
uncertainty in Alog g yields excessively large error bars on 
our determination of the radius and, consequently, the mass, 
and for this reason we deem this method unreliable with our 
current data set. A second method proves to be more fruitful 
and is based upon the Baade-Wesselink method. Any frac- 
tional change in the angular radius should yield the same 
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Figure 11. Plot of Aa vs AR. For a radial pulsator, these vari- 
ables should exhibit a linear relationship. The dashed line denotes 
the linear fit to these data. The data exhibit a Pearson correla- 
tion coefficient of 0.78. A representative error bar is shown for 
the leftmost data point and was computed by propagating errors 
through Eqn 3. 



fractional change in the physical radius. Using this fact and 
linearizing Eqn. 3, we compute R from the following relation: 



AR 

R 



Aa 



1 

2 V h 



AF\ 

Fx 



(5) 



0.19 ± 0.08 Rq. Sim- 



This expression yields a radius of R : 
ilar results are obtained by calculating the radius directly 
from a and d, which gives R = 0.20 ±o:o6 ^0- We present 
the derivation of R using Eqn 5 to illustrate that errors in the 
observed ffux cancel out and do not affect the computation 
of the physical radius as they do the distance calculation. 

Combined with our measurement of the mean surface 
gravity, the radius derived from Eqn 5 gives a mass of M 
= 0.39 ±o;i3 Mq. Again, the radius and mass we report 
are meaningless if the pulsation is non-radial. Although our 
errors bars are large, the mass agrees with the commonly- 
accepted value of 0.5 Mq. The radius is also consistent with 
those determined for other sdB stars from binary obser- 
vations and from theoretical mass-radius relationships but 
has been computed using more direct observations. Work- 
ing backwards and assuming a canonical mass of 0.5 Mq, 
we would derive a radius of R = 0.22 Rq from the mean log 
g value. 



5 CONCLUSIONS 

We report the discovery of an exciting new member of the 
sdBVr class of pulsators, CS 1246. Our simultaneous, multi- 
colour, time-series photometry reveals a single oscillation 
mode with period 371.7 s. The amplitude of this mode de- 
pends strongly on wavelength and is largest in the u' filter. 
We detect no additional periodicities in the frequency spec- 
tra, although we note our 0.41-m telescope detection limits 
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are poor. The light curve of CS 1246 appears to be the sim- 
plest of all sdBV stars. 

Time-series spectroscopy proved to be a much more 
fruitful exercise than we originally anticipated. We detected 
variations in temperature, gravity, and radial velocity at 
the same frequency as the photometric oscillations. Inter- 
estingly, the ratios of these amplitudes (507 K, 0.034, 8.8 
km s~^) scale extremely well with those found by 0stensen 
et al. (2007) for the main mode of Balloon 090100001 (1186 
K, 0.084, 18.9 km s"^), which Baran et al. (2008) have con- 
firmed as a radial oscillation. By comparing changes in the 
angular and physical radii, we have shown the pulsation 
mode in CS 1246 is consistent with a radial one, but we 
cannot rule out non-radial modes with our current data set. 
We were able, however, to calculate the radius and mass 
under the assumption of a radial pulsation, and find values 
consistent with those derived for other sdB stars. Our study 
is the first to apply a method based on the Baado-Wesselink 
technique to an sdBV star. There is room for improvement 
from long time-series spectroscopy runs, though, since the 
uncertainty in our mass does not allow one to distinguish 
between various sdB formation scenarios. Better photome- 
try and more precise radial velocities arc required to achieve 
such a goal. If the 371. 7-s mode is confirmed to be radial, 
application of the Baade-Wesselink method to a more sub- 
stantial data set could provide the most direct measurements 
of an sdB radius and mass ever made (except in cases where 
the sdBV is in an eclipsing binary system). 

The atmospheric parameters we derived place CS 1246 
near the boundary of the sdBVr and sdBVs instability 
strips, and, consequently, it could be a hybrid pulsator. To 
date, three such sdBVrs stars have been noted in the liter- 
ature: Balloon 090100001 (Orciro ct al. 2005; Baran et al. 
2005, 2006), HS 0702-1-6043 (Schuh et al. 2005, 2006; Lutz 
et al. 2008), and HS 2201-1-2610 (Lutz et al. 2009). Since 
studies of rapid and slow oscillations probe different depths 
in a star, the sdBVrs stars can provide much more informa- 
tion than pulsators exhibiting only one of these modes. The 
amplitudes of slow oscillation modes tend to be small (a few 
mma), and we cannot rule out the possibility that CS 1246 
is a hybrid pulsator since our detection limits in the low- 
frequency regime (2-3 mma) are comparable to this level. 
A longer photometry run with a larger-aperture telescope is 
needed to search for slow pulsation modes. 

The study of CS 1246 is complicated by the presence 
of the Coalsack Dark Nebula along our line of sight, which 
helps to explain two anomalous observations. The first prob- 
lem is the B-V colour reported by Reid et al. (1988), which 
is too red for an sdB star. Previous extinction studies of the 
Coalsack report reddening values similar to the one we find 
for CS 1246, so it is not impossible interstellar extinction 
is the only source of the observed reddening. The Coalsack 
might also explain why our derived distance is less than 
what one would expect given the brightness of the star. The 
effective temperature and radius we calculated imply an ab- 
solute Visual magnitude of 4.4 mag. Combining this absolute 
magnitude with the apparent one, we derive a distance of 
nearly 1100 pc, much farther than the 460 pc calculated in 
§4.7. Accounting for the approximately 1.8 mags of Visual 
extinction expected from the Coalsack, however, makes our 
computed distance consistent with the apparent magnitude. 
We also note that at 460 pc, CS 1246 is located well behind 



the nebula, which has a distance of 174 pc (Rodgers 1960). 
If not for the Coalsack Dark Nebula, the star would be one 
of the brightest known sdBVs. 

Finally, we draw attention to the potential CS 1246 has 
for the detection of its secular evolution and the discovery 
of orbiting low-mass companions. Since the frequency spec- 
trum is dominated by a single, large-amplitude oscillation, 
phase shifts should be relatively simple to detect through an 
0-C diagram. The shifts may come about from secular evo- 
lution of the star as it evolves along or off of the extended 
horizontal branch or from wobbles due to orbital interac- 
tions with unseen, low-mass companions. This method has 
already been used successfully to find a planet around the 
sdB star V391 Pcgasi (Silvotti ct al. 2007). As PROMPT 
is an ideal tool for such a study, we are currently building 
up an ephemeris for CS 1246 and look forward to producing 
our first 0-C diagram in the future. 
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APPENDIX A: CORRELATION COEFFICIENT 
PROBABILITIES 

We compute the correlation coefficient and corresponding 
probabilities using the methods described by Taylor (1997). 
The linear correlation coefficient r between N points (xi,yi) 
can be computed using the expression 



(Al) 



The coefficient r will be close to ± 1 if the linear correlation 
is strong and near if the points axe uncorrelated. As even 
uncorrelated points can produce a non-zero coefficient, one 

must further quantify the significance of r by considering the 
number of measurements. The probability that N measure- 
ments of two uncorrelated variables result in a coefficient 
greater than or equal to ro is given by 



PN{\r\ > \ro\) 



2r[{N - l)/2] 
^r[{N - 2)/2] 



(l-r^) 

Jlr„l 



,2n(JV-4)/2 



dr.(A2) 



This probability strongly depends upon the number of mea- 
surements N; the probability that two uncorrelated variables 
give a set correlation coefficient decreases with increasing N. 

In our case, the Aa and AR points measured give a 
coefficient of 0.78. The above expression shows there is less 
than a 2% chance that two uncorrelated variables give a 
cocflicicnt at or above 0.78. Thus, wc have significant but 
not definitive evidence of a linear correlation. 

This paper has been typeset from a T^JC/ I^TJ^X file prepared 
by the author. 



